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Abstract 

The Large Synoptic Survey Telescope (LSST) will use an active optics system (AOS) to maintain alignment 
and surface figure on its three large mirrors. Corrective actions fed to the LSST AOS are determined from 
information derived from 4 curvature wavefront sensors located at the corners of the focal plane. Each 
wavefront sensor is a split detector such that the halves are 1mm on either side of focus. In this paper we 
describe the extensions to published curvature wavefront sensing algorithms needed to address challenges 
presented by the LSST, namely the large central obscuration, the fast //1.23 beam, off-axis pupil distortions, 
and vignetting at the sensor locations. We also describe corrections needed for the split sensors and the 
effects from the angular separation of different stars providing the intra- and extra-focal images. Lastly, 
we present simulations that demonstrate convergence, linearity, and negligible noise when compared to 
atmospheric effects when the algorithm extensions are applied to the LSST optical system. The algorithm 
extensions reported here are generic and can easily be adapted to other wide-held optical systems including 
similar telescopes with large central obscuration and off-axis curvature sensing. 
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1. Introduction 


The Large Synoptic Survey Telescope (LSST) is a new facility now under construction that 
will survey ~20000 square degrees of the southern sky through 6 spectral Liters [ugrizy] 
multiple times over a 10-year period [T]. The conduct of the LSST survey is dehned by a 

2- exposure “visit” lasting ~39 seconds having the following pattern: 15s integration plus 
Is shutter open/close transitions (16s elapsed time), 2s focal plane array (FPA) readout, 
16s exposure, 5s telescope repointing and FPA readout (the last 2s readout is concurrent 
with the telescope repointing). With this cadence the LSST will observe ~800 visits per 8 
hour night, producing more than 1600 images and ~15 Tbytes per night of data and cover 
the entire visible sky every 3-4 nights. Construction of the LSST includes an 8-meter class 

3- mirror telescope, a 3.2 billion pixel camera, and an extensive computing system for data 
analysis and archiving. Commissioning of the LSST is expected to begin in late 2019 and 
full surveying operations at the end of 2022. 

The optical system (Figure left) of the LSST is based on a modihed Paul-Baker 3- 
mirror telescope design having an 8.4m primary, 3.4m secondary, and 5.0m tertiary feeding 
a three-element refractive camera system producing a flat 3.5 degree field-of-view with an 
effective clear aperture of 6.5m. The mean intrinsic imaging performance of the optical 
design across the optical spectrum from 320nm to 1150nm is better than 0.1 arcsecond full 
width at half maximum (FWHM) over the full held of view. The optical design has allowed 
the primary and tertiary mirrors to be fabricated out of a single substrate. Fabrication 
of the primary-tertiary (M1M3) mirror pair has recently been completed by the Steward 
Observatory Mirror Lab. 

The overall system image quality budget for the LSST is 0.4 arcsecond FWHM and is 
allocated between the telescope (0.25”) and camera (0.3”) subsystems. With this image 
quality budget the LSST’s delivered image quality is dominated by the atmospheric seeing 
at its site on Cerro Pachon in Chile. For the telescope subsystem, the bulk of its budget 
allocation is taken up by residual mirror hgure errors and thermal and gravity-induced 
misalignments of the camera and secondary mirror systems. In order to maintain this 
image quality budget the LSST telescope subsystem utilizes an active optics system (AOS) 
controlling 50 degrees of freedom (DoF) consisting of 20 bending modes each on the actively 
supported M1M3 and M2 systems and 5 DoF of rigid body position for the camera and M2. 
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Fig. 1. The 3-mirror telescope optical configuration for the LSST (left) showing the placement 
of the tertiary mirror within the primary mirror, allowing the two to be fabricated from a single 
substrate. The camera optics (right) use three fused-silica lenses and a meniscus filter. The r-band 
image quality (lower right) of the LSST optical system is shown for 0, 1.0 and 1.75 degrees from 
the optical axis. 

By contrast, in the camera snbsystem mnch of its bndget allocation is taken up by internal 
charge diffusion of the 100/rm-thick CCDs and overall flatness of the FPA. 

Most modern astronomical telescopes use some form of AOS with near-realtime optical 
wavefront sensing feedback. The most common method for estimating the wavefront in 
these systems is with a Shack-Hartmann wavefront sensor (SHWFS), e.g. the Magellan 
telescopes [ 2 ] and the VLT [ 3 ]. Wide field systems with their fast optical beams make using 
Shack-Hartmann wavefront sensing problematic for two reasons: 1) To avoid the vignetting 
effects a pickoff mirror would cause in a fast converging optical beam, the SHWFS would 
have to be placed near the focal surface. This means these mechanisms would have to be 
small and within the vacuum cryostat of the imager. 2) Some form of articulation is needed 
to place the reference star within the SHWFS limited field of view. 

Currently operating wide held survey telescopes all use wavefront sensing based on image 
analysis from area detectors. The VISTA telescope, an infrared survey telescope, utilizes a 
pair of curvature sensors and forward modeling of low order aberrations to provide optical 
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feedback W- Similarly, the Dark Energy Camera employs 8 2k x 2k sensors around the focal 
plane both in and out of focus to forward model the resulting donut images [U]. Curvature 
wavefront sensing is also in use on several operating telescopes including the 3.5m WIYN [7] 
and 4m Mayall [8] telescopes at Kitt Peak National Observatory. These telescopes use 
curvature wavefront sensing to build their AOS look-up tables and to periodically set the 
operational zero points for the AOS. In each of these systems the number of DoF controlled 
is limited to focus, misalignment of the instrument or secondary mirror (coma), and one or 
two forms of primary mirror bending (astigmatism and/or trefoil). 

Because LSST uses two actively supported mirror systems and two positioning hexapods, 
the number of controlled DoF is signihcantly greater than typically seen in currently operat¬ 
ing telescopes. Controlling these DoF requires the ability to estimate higher order properties 
of the aberrated wavefront, Zernike coefficients Z4 - Z22 in Noll/Mahajan’s dehnition [^1X0]. 
Further, the time constraints imposed by the LSST’s rapid cadence requires that the process 
of measuring the optical wavefront and maintaining alignment and hgure control be highly 
automated and reliable. For these reasons we have adopted curvature wavefront sensing to 
provide the optical feedback needed to control the LSST’s AOS. 

The layout of the LSST focal plane (Figure]^ left) was developed in part to accommodate 
the LSST’s wavefront sensing requirements. The science sensors (blue) are arranged in 21 
modular “rafts”, each containing 9 CCDs. At four locations are “corner rafts” consisting 
of two guide sensors (yellow) and one wavefront sensor (green) each. By measuring the 
held-dependent Z4 - Z22 at 4 locations we have 76 variables to control 50 DoF. 

Curvature wavefront sensing offers some advantages over other wavefront sensing methods 
for wide held survey telescopes. By relying on equally defocused intra- and extra-focal 
images, curvature sensing can use area sensors with relatively large £elds-of-view. This allows 
significant hexibility in selecting reference sources to use for the wavefront measurement 
when the source scenery is constantly changing from one visit to the next with the LSST 
cadence. 

Typically this is achieved through using a beam splitter and delay line or by physically 
moving the detector. Both of these approaches are problematic for the LSST with its fast /- 
number (//1.23) and crowded focal plane. Our design for LSST therefore splits the available 
wavefront sensor area in two halves (see Figure]^ right), with one 1mm in front of nominal 
focus and the other 1mm behind. Each half utilizes a 2k x 4k CCD with lO/im square pixels 


4 


Wavefront Sensors 




CCS ♦X 



Fig. 2. The focal plane configuration for the LSST (left panel) showing sensors used for science 
(blue), guiding (yellow), and wavefront sensing (green). Each raft is mounted into a preadjusted 
grid “bay”. The wavefront sensors are divided into two halves of intra- and extra-focal sensors 
(right panel). The strategy in making use of multiple stars on each half-chip to obtain a wavefront 
measurement at each corner is currently under investigation and not discussed in this paper. 

providing a 7x14 arcminute field of view. With this field of view, the likelihood of acquiring 
suitable reference stars is near unity even in the low density galactic pole regions m- 
Therefore, no active acquisition of the reference source is required. 

There are a variety of algorithms used to estimate the wavefront in curvature wavefront 
sensing. However, none of them by themselves can work with a wide-held telescope having 
large central obscuration and off-axis sensors like LSST. Our strategy is to choose two 
well-established algorithms that are known to work for large /-number, on-axis systems, 
implement them as our baseline, then extend them to work with small /-number and off- 
axis sensors. As our two baseline algorithms we have chosen the iterative Fast Fourier 
Transform (FFT) method by Roddier and Roddier [ 12 ], and the series expansion technique 
by Gureyev and Nugent na. We found both methods to be accurate and reasonably fast. 
We note that our extensions to the baseline algorithms can also be used with other curvature 
sensing algorithms that work for large /-number and on-axis systems, in order to make them 
work for small /-number and off-axis sensors. 
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In this paper, we describe the development of the curvature wavefront sensing algorithm 
that measures the wavefront at the 4 corner locations. The paper is organized as follows. In 
Section we review our baseline algorithms of curvature wavefront sensing, and in Section 
we discuss the new challenges facing a wide-held and off-axis system like LSST. We also 
discuss in Sectionj^the required modihcations to the baseline algorithms in order to overcome 
these challenges. Section then gives simulation results of unit testing and validations, 
including analyses of algorithmic noise and atmospheric background. 

2. Baseline Curvature Sensing Algorithms 

The concept of curvature wavefront sensing was hrst developed and demonstrated by F. 
Roddier [TT]. The underlying idea is to measure the spatial intensity distribution of a star 
at two positions, one on either side of focus. The derivative of the local surface brightness of 
the defocused images along the direction of propagation is given by the transport of intensity 
equation (TIE) [12]: 

dl 

— =-(VJ-VIT + JV^IT), (1) 

where I is the intensity, W is the wavefront error, and 2 ; is the distance between the conjugate 
planes of the defocused pupil images. 

Assuming the intensity is equal to Jq everywhere inside the pupil and zero outside, we 
have the Neumann boundary conditions: 

VI = -Ioh6,, (2) 

where 6c is a delta function around the pupil edge. Therefore, the hrst term in Eq. ([^ is 
localized at the beam edge. Across the beam, the derivative of the local image brightness is 
proportional to the Laplacian, or curvature, of the wavefront. 

As a partial linear differential equation, the TIE can be solved with various methods. In 
the following 2 subsections we summarize the two methods we have adopted as our baseline 
algorithms. 

2.A. Iterative FFT 

The iterative FFT algorithm na solves the TIE by making use of the Laplacian operator 
becoming a simple arithmetic operation in Fourier space. 
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The longitudinal derivative of the intensity can be expressed as 


1 dl 
lodz 


= - 5. 


dW 

dn 


( 3 ) 


We define the wavefront signal S using the longitudinal derivative normalized by Jq. As 
such, S can be approximated as: 


S = -i^ 

lodz 


h-h 


1 h-h 


h ■ Az h T h 

where the focus offset in the object space (e.g. the pupil) is given by 


( 4 ) 


Az = f{f-l)/l, 


( 5 ) 


with / the system focal length, and I the defocus distance of the intra/extra focal image 
planes. 

It can be shown that if we constrain dW/dn on the edge, 6c can be absorbed into the 
Laplacian. An estimate of the Laplacian of the wavefront error can be rewritten as 

^ S. (6) 


Using the following property of the Fourier Transform (FT) (see for example page 314 of 
Ref. [H]) , 

FT{fi,u){V^W{x,y)} = + 

FT{fx,u){W{x,y)}, (7) 


where y and iz are the spatial frequencies, we can solve for W using the inverse Fourier 
Transform {IFT): 


W K, IFT{x,y) 


FT{y,iy){S} 


( 8 ) 


-47r2(p2 + 

The implementation of the FFT algorithm is shown in Figure (left, boxes 1-8). It 
involves iteratively applying Equation |^to estimate the wavefront W, then setting boundary 
condition and putting the original wavefront signal, S, back inside the boundary to obtain 
a new signal. The algorithm we have described thus far is called the “inner loop”. 

As was noted by Roddier and Roddier ns, the method above is only a first order approx¬ 
imation valid for small values, i.e., highly defocused images (see Eq. ([^). The initial 
solution of the Poisson equation we obtain by the “inner loop” should be used as a first 
order solution that is further rehned in a second iterative process. 
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Fig. 3. The block diagram of the iterative FFT algorithm including the “outer loop” image 
compensation for the estimated wavefront aberrations. 

The overall algorithm accuracy can be improved by iteratively removing or compensating 
for the effects of the estimated wavefront aberrations on the original intra- and extra-focal 
intensity images and then reapplying the “inner loop” to the corrected images. This “outer 
loop” (Figure (right, boxes a-e)) is iterated on until the noise level or a given number 
of iterations is reached. With each iteration of the “outer loop” the estimated residual 
wavefront is summed with the previous solution to provide the current wavefront estimate. 

We compensate/correct the intra- and extra-focal images for the aberrations of the current 
wavefront estimate by remapping the image flux using the Jacobian of the wavefront in the 
pupil plane. For this process we let R be the pupil radius of the telescope and denote the 
reduced coordinates in the pupil plane as a; = U/R and y = V/R, where U and V are 
Cartesian coordinates in the pupil plane, and those in the image plane as x' and y'. The 
reduced coordinates in the pupil and image planes are related by na 


x' = X + CdW {x, y) / dx^ 
v' = y + CdW(x,y)ld-ii, 


( 9 ) 

( 10 ) 
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( 11 ) 


c = 


The intensities in the two planes are related by the Jacobian (J) due to flux conservation, 


dx'/dx dx'/dy 
dy'/dx dy'/dy 

Using these relations, given a wavefront estimation from a pair of defocused images, we 
are able to “restore” the pair of images to a state where the given estimated wavefront 
aberrations are absent. The new images are then fed to the inner loop (see box c in Fig¬ 
ure where the residual wavefront error is estimated. The loop continues until the residual 
wavefront JFresiduai is consistent with zero, or a given number of outer iterations is reached. 

Our compensation algorithm has an oversampling parameter, which enables sub-pixel 
resolution for the mapping between the pupil and image planes. We can choose to improve 
the compensation performance by increasing the sub-pixel sampling, but at a cost of com¬ 
putation time. Due to the fact that the compensation is based on geometrical ray-tracing, 
we always compensate on the original defocused images. A feedback gain less than unity 
is used to prevent large oscillation in the final wavefront estimation, i.e., upon each outer 
iteration only part of the residual is compensated. To decouple Zernike terms with the 
same azimuthal frequency, for example, between tip-tilt and astigmatisms, and defocus and 
spherical aberration, only a certain number of low order Zernike terms are compensated at 
each outer iteration. After the algorithm has converged on low order terms, higher order 
terms are added to the compensated wavefront. The tests we show in Section |^each include 
14 outer iterations, with the highest Zernike index of 4, 4, 6, 6, 13, 13, 13, 13, 22, 22, 22, 
22, 22, and 22, respectively. 

The compensation algorithm also helps identify cases where the geometric limit has been 
reached. When the compensation procedure results in a negative intensity on the pupil 
plane, an in-caustic warning flag is set. Here “in-caustic” means rays approaching focus 
from different points in the pupil cross before the intra focal or beyond the extra focal image 
plane. This crossing leads to ambiguities in the interpretation of the image intensities and 
can lead to non-physical negative intensities in the compensated image. 



Hx,y)/I'{x,y) 
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2.B. Series Expansion 


The series expansion method of solving the TIE is based on the decomposition of the TIE 
into a series of orthonormal and complete basis fnnctions [ 13 ] . Since the onter loop is nsing 
the annnlar Zernike polynomials, it is natnral to choose those for the basis, albeit the method 
wonld certainly work with any orthonormal basis set. 

Let n be the area of the plane with positive intensity I and smooth bonndary T. As 
snch, the TIE can also be written as: 


dj = -V ■ (/VhE) (13) 

I{x,y) > 0 inside hi (14) 

l{x,y) = 0 ontside n and on T. (15) 


Let Zi{x,y), i=l,2,3... be a set of orthonormal and complete basis fnnctions over the pnpil, 

OO 

W{x,y) = J2WiZ^{x,y)■ (16) 


2=1 


Now we mnltiply Eq. (13) by Zj{x,y) and integrate over hi. 




{dj)Zjdn = - / V ■ {ivw)Zjdn. 


( 17 ) 


Integrating by parts, and by taking into acconnt the boundary condition, we get 


where 


Therefore, 


MjiWi = Fj, 

2=1 





{dj)Zjdn 


Mji 


I{x,y)'VZj ■ VZidVL. 


W = M^^F. 


(18) 


(19) 

( 20 ) 


( 21 ) 


Each of our baseline algorithms have certain advantages and applications over the other. 
The wavefront compensation we have discussed for the iterative EFT method can also be 
used to improve the performance of the series expansion technique, or any other algorithm 
used to solve the TIE. For detailed analysis of the wavefront structure the iterative EFT is 


10 



preferred since it results in a 2-D image of the wavefront. Both methods work with arbitrary 
pupil geometry, provided that a set of orthonormal basis functions over the pupil can be 
found. When such basis functions are not available, the accuracy of the series expansion 
method degrades more because the expansion relies on the orthogonality of the basis, whereas 
for iterative FFT, the non-orthogonality only makes decomposition of the solved wavefront 
into the Zernike space problematic. On the other hand, when speed is a concern, the series 
expansion is the preferred method. With the same number of outer iterations, we found 
that the series expansion is about 5 times faster than the iterative FFT. We anticipate both 
methods will be used in LSST during commissioning and engineering (iterative FFT) and 
routine survey operations (series expansion). In the analyses presented in the rest of the 
paper, we use the series expansion as the inner-loop, and the wavefront compensation as the 
outer-loop. 


3. Algorithm Modifications for LSST 

LSST’s optical system poses four algorithmic challenges to using curvature sensing and 
solving the intensity transport equation for estimating the wavefront error at the entrance 
pupil: 1) the high central obscuration requires the use of annular Zernike polynomials as 
the basis set; 2) the fast //1.23 optical beam results in signihcant nonlinearity in projecting 
the wavefront error on the pupil plane; 3) at the location of the wavefront sensors there is 
signihcant pupil distortion and vignetting that must be accounted for; and 4) held-dependent 
variations in the wavefront over the wavefront sensor area must also be accounted for. Our 
solutions for each of these are given in the discussion that follows. 


3.A. Large Central Obscuration 

In both the iterative FFT and the series expansion methods, we decompose the wavefront 
onto a set of orthonormal and complete basis functions. Traditionally, the “standard” hlled 
aperture Zernike polynomials have been chosen for this purpose due to their well-understood 
properties [9]. When the optical system has non-zero obscuration, the Zernike polynomials 
are no longer orthogonal. When the obscuration is small, the non-orthogonality does not 
cause much of a problem. However, when the central obscuration gets large, the non- 
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Fig. 4. Orthogonality test of the standard Zernike polynomial on an entrance pupil with 60% 
obscuration. The x- and y- axes are the order of the Zernike polynomials (up to 22). 

orthogonality causes a degeneracy problem in the Zernike decomposition. Figure shows 
the orthogonality test of the standard Zernike polynomial on an LSST-like 60% obscured 
entrance pupil. 

The obvious solution is to use the annular Zernike polynomials [10]. This ensures the 
orthogonality on the annular pupil. However, note that an obscuration means there is a 
loss of information. Even with the annular Zernike polynomials, large central obscuration 
means that the wavefront is harder to measure compared to the case when there is no or 
small obscuration. 

3.B. Correcting for Fast /-number 

In a fast-beam optical system like LSST the wavefront error dehned on the reference sphere 
at the exit pupil can no longer be projected on the pupil plane in a straightforward, linear 
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way. The steepness of the fast-beam reference sphere results in non-linear mapping from the 
pupil {x — y plane) to the intra- and extra-focal images (x' — y' plane). Consequently, the 
intensity distribution on the aberration-free image is no longer uniform and would give rise 
to anomalous wavefront aberration estimates. This effect also causes the obscuration ratio 
as seen at the intra- and extra-focal image planes to be different from that on the pupil; 
thus the linear equations for the mapping between x — y and x' — y' (Eqs. and 10) are 
no longer valid. Therefore, we need to quantify the non-linear fast-beam effect, and be able 
to remove it so that the corrected images can be processed using the iterative FFT or the 
series expansion algorithms. 

For an on-axis fast-beam system, it can be shown that 


x' = F(x, y)x + C 
y' = F{x,y)y + C 


dW 

dx ’ 
dW 


( 22 ) 

(23) 


with 


F{x,y) = 


m 




\Jp - {x‘^ F y‘^)F?'' 


(24) 


where m = R' f /{IR) is the mask scaling factor, R' is the radius of the no-aberration image, C 
has been dehned in Eq. (0, and the reduced coordinates on the image plane are normalized 
using the paraxial image radius R!/m = IR/f, i.e., the radius of the image for a paraxial 


system with the same / and R. Note that flux conservation (Eq. (12)) still holds. The 
partial derivatives of the Jacobian become 


dx' 

dx 

dy' 

dy 

dx' 

dy 

dy' 

dx 


r./ Q2^r 

+ p_^.2^2)+^ Q^2 

(25) 

Pr^ d‘^W 

(26) 

f2 _ ^2^2 + ^ Q^y 

(27) 

y^ f 2 _ r 2 R 2 dxy' 

(28) 


With this mapping we are able to fully account for geometric projection effects of the LSST’s 
fast /-number. 
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3.C. OfF-axis distortion and Vignetting Correction 


The LSST wavefront sensors are located about 1.7° off the optical axis. On top of the 
non-linear projection effects due to the small /-number, the off-axis images are distorted 
and non-axisymmetric. Since an analytical mapping from the telescope aperture to the 
defocused image planes is not readily achievable, we instead have developed a numerical 
solution, representing the mapping between the two sets of coordinates with 2-dimensional 
10th order polynomials. The coefficients of the polynomials are determined using least- 
square fits to the ray-hit coordinates from a grid of rays simulated by the LSST ZEMAX 
model. The agreement between the fit and the ZEMAX data is better than one percent of 
a pixel. The gradients and Jacobians are then calculated and implemented in the wavefront 
compensation code. 

Compared to a paraxial model with the same /-number, the off-axis distortion shifts the 
ray-hit coordinates on the image planes by up to about 6 pixels, relative to the chief ray, 
had there been no central obscuration. At the upper right corner of the focal plane, this 
distortion is symmetric about the 45° line. In contrast, the maximum shifts in the ray-hit 
coordinates on the image plane caused by 200nm of 45°— astigmatism is about 0.2 pixel. 

Because the wavefront sensors are at the edge of the 3.5° field of view, the wavefront 
images are vignetted. For most parts of the wavefront chip, vignetting is mainly due to 
the primary and secondary mirrors and increases gradually with field angle. When the field 
angle gets to about 1.7°, vignetting due to the camera body starts to cause a sharp decrease 
in the fraction of unvignetted rays. Figure j^shows intra-focal images on a 5x5 grid covering 
the upper right wavefront chip with field angle ranging from 1.51° to 1.84°. 

Vignetting means loss of information on the edge of the pupil, making it harder to recover 
the wavefront. This is especially true for aberrations like coma, where most of the intensity 
variations occur near the edge of the image. Furthermore, although neither the iterative 
FFT nor the series expansion algorithm has requirements on the shape of the pupil, they do 
rely on decomposition of the wavefront onto the annular Zernike polynomials, which are no 
longer orthogonal on vignetted pupils. Because the vignetting is relatively small, ~10%, the 
deterioration of the wavefront estimation accuracy due to this non-orthogonality has been 
observed to be small. At the center of the wavefront chip, unit tests using ZEMAX images 
show no visible increase in estimation uncertainty compared to the on-axis tests. 
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Fig. 5. Intra-focal images on a 5x5 grid covering the LSST upper right wavefront chip with field 
angle ranging from 1.51° to 1.84°. Note that these are ZEMAX simulated images for analysis and 
testing purposes. In reality, the intra-focal images are acquired on half of the wavefront chip only. 

3.D. Field-Dependent Corrections 

The split sensor design of LSST introduces two additional potential sources for error. First, 
the TIE used in curvature sensing requires that the total flux in the intra- and extra-focal 
images are the same. Using different sources for the intra- and extra-focal images in general 
will violate this requirement. This can be largely overcome by local background subtraction 
and normalizing to the source flux during the processing of the intra- and extra-focal images 
prior to solving for the wavefront. 

Second, the curvature sensing method assumes that both the intra- and extra-focal images 
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contain the same optical properties (aberrations and pupil geometry). Because the intra- 
and extra-focal sources will come from different parts of the held of view they will not have 
identical optical properties. 

The off-axis geometric distortions are held-dependent as are the telescope intrinsic aber¬ 
rations and vignetting. The off-axis correction we discussed previously corrects intensity 
variations, regardless of their source, be it from the off-axis distortion or the held-dependent 
telescope intrinsic aberrations. When we make the oh-axis correction, we always do so ref¬ 
erenced to a common held coordinate at the center of the wavefront sensing area, regardless 
of the source’s original location. Thus, the oh-axis correction removes both the ehects of 
held dependent geometric distortion and intrinsic design aberrations. 

The vignetting pattern still varies by held. To avoid anomalous signal due to the diherence 
in vignetting from entering the wavefront estimation results, we mask oh the intra- and extra- 
focal images using a common pupil mask, which is the logical AND of the pupil masks at 
the intra- and extra-focal held positions. 


4. Algorithm Validation and Unit Testing 

We have applied many levels of unit testing on the LSST wavefront sensing algorithms using 
images with and without the ehects of atmospheric turbulence. Examples of these tests and 
their results are provided below. 


4.A. Paraxial Unit Testing 

We have performed comprehensive tests with on-axis images created in ZEMAX using a 
paraxial lens model, with /-number ranging from 1.3 to 4; image space defocus from 1mm 
to 5mm; and obscuration ratios from 0 (hlled) to 60%. One example of such tests is given 
in Figure]^ The results agree nicely with the ground truth from ZEMAX. 

Our wavefront estimation algorithm performs very well. For the single-Zernike-term 
tests, each Zernike mode runs into caustic at a diherent aberration magnitude. For //1.23, 
1mm image space defocus, with RMS wavefront aberration below 1.5 microns, the only 
problematic modes are the 5th order astigmatisms. 
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4.B. LSST On- and Off-axis 


Unit test results using ZEMAX simulated images for various Zernike terms at different 
magnitudes show that the fast-beam correction algorithm is very successful. To facilitate 
algorithm testing, the telescope intrinsic aberrations have been compensated using a phase 
screen, before applying the known aberration under test. 

One example of such tests is given in Figure The intra- and extra-focal images are 
taken at the center of the focal plane. In this example, the aberration present is 0.5 wave 
(wavelength is 770nm) of Zll (spherical aberration). The estimated Zll, as shown in Fig¬ 
ure agrees nicely with the ground truth. All other Zernike terms have magnitudes of no 
more than 20-30nm. 

We have also performed unit tests using ZEMAX simulated LSST off-axis images with 
various Zernike terms at different magnitudes. One example of such tests, where 2 waves 
(wavelength is 770nm) of Z5 (45° astigmatism) is applied, is shown in Figure]^ The intra- 
and extra-focal images are taken at held position (1.185°, 1.185°), corresponding to the 
center of one of the LSST wavefront sensors. Vignetting is clearly visible in these images. 
Again, the input aberration is recovered, as shown in Figure All other Zernike terms 
have magnitudes very close to zero. Tests are also performed using images simulated by 
the LSST Photon Monte Carlo Simulator im. A similar level of agreement between the 
measured wavefront and the truth is observed. 

As another set of tests, the ZEMAX model is used including all intrinsic aberrations and 
various individual perturbations at the 4 corner wavefront sensor locations. Figure shows 
an example at the held angle of (—1.185°, 1.185°), where the secondary mirror has been 
decentered by 0.5mm. 


4.C. Algorithm Linearity 

By varying the controlled degrees of freedom of the telescope model and repeating the 
wavefront measurements, good algorithmic linearity has been observed within the geometric 
limit. One example of such tests is shown in Figure When the secondary mirror is tilted 
by incremental amounts, and all other degrees of freedom of the telescope stay unperturbed, 
the wavefront Zernikes (45°-astigmatism shown) change linearly. The linearity starts to 
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degrade in proximity of the geometric limit. Similar behavior has also been observed for 
other degrees of freedom of the system. 


4.D. Covariance Analysis 


The atmospheric covariance is computed for 15-second time integrated atmospheres gen¬ 
erated using the Arroyo library na. Six layers of Kolmogorov phase screens at various 
heights are simulated. Based on historic DIMM data at Cerro Pachon, and the current 
understanding on its outer scale, we use tq = 17cm and assume inhnite outer scale. 

The same time-integrated phase screens used for the atmosphere covariance are then 
applied to a ZEMAX ray trace model in order to evaluate the importance of algorithm noise 
with respect to atmosphere noise in estimating the wavefront. For each time-integrated 
phase screen instance, ZEMAX is used to generate intra- and extra-focal image pairs at 
each of the 4 wavefront sensors. These images are then processed through the LSST specihc 
curvature algorithms to estimate the wavefront coefficients. Compared to the calculated 
coefficients from atmosphere alone (ideal sensor), the results from estimating the wavefront 
using the simulated sensors are similar, indicating that the covariance is dominated by the 
atmospheric contribution and not the wavefront sensing algorithms. The plots in Figure 10 
show the diagonal entries (19 entries for each held point, corresponding to the Zernike terms 
estimated) and the singular values of the covariance matrix; the largest singular values most 
important for estimating the misalignment of the telescope are almost identical. 


5. Summary, Conclusions, and Future Work 

Extensions to two curvature wavefront sensing algorithms for LSST have been developed. 
The underlying algorithms considered are the iterative EFT method by the Roddiers [12] and 
the series expansion based method by Gureyev and Nugent [T3|. These are well-established 
methods that have been proven to work well for paraxial systems. Several modihcations are 
needed to make them work for LSST, to overcome challenges including a highly obscured 
pupil, the fast /-number, pupil distortion and vignetting at the held corners, and variation 
of the wavefront over the area covered by the split-sensors. Our baseline algorithm for use in 
routine operations is the series expansion method, due to its higher computational efficiency. 
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Extensive simulations have been performed using images generated by ZEMAX and the 
LSST Photon Monte Carlo simulator [T7] . 

Integrated modeling of the LSST AOS, including the optimal control of the system, as 
well as adapting the algorithms to be used by several other operating telescopes, is reported 
in [in]. As the LSST construction progresses we are also looking into alternative wavefront 
sensing algorithms, such as the forward modeling technique used by DECam, and the PSF- 
based technique used by VST, to compare their performance to the results shown here. 

This is the hrst of a series of papers being planned on the Active Optics System of LSST. 
The following topics are planned to be discussed in subsequent papers hence not included 
here: 1) the validation of the algorithms presented in this paper using real data taken 
at major operating telescopes; 2) the LSST alignment strategy, the tomographic optical 
reconstruction, and the telescope control algorithm; 3) the general functionality of the LSST 
active optics systems (AOS) and the prototype processing pipeline we are using to develop 
a full simulation of the LSST AOS operation. 
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Fig. 6. An example of a paraxial lens unit test for focus (Z4). The top row shows the input 
images. The middle row shows how the estimated Z4 converges to the input value. The final 
annular Zernike composition is shown at the bottom. Note that on the bottom plot, we show 
results from the last hve outer iterations. These curves are on top of each other, which shows good 
convergence. 
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Fig. 8. ZEMAX simulation at the field angle of (—1.185°, 1.185°) containing all intrinsic aberra¬ 
tions plus a decenter of 0.5mm in the secondary mirror. The recovered aberration coefficients (red 
circles) are in excellent agreement with the ZEMAX-calculated coefficients (blue crosses). 
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Fig. 9. ZEMAX-calculated wavefront (blue crosses) and the recovered wavefront (red circles) at 
one of the four corners of the focal plane. The recovered aberration coefficients start to deviate 
from linearity in proximity of the geometric limit. 
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Fig. 10. The nearly identical covariance and singular values of the ideal sensor (atmosphere only) 
and fully simulated and processed images show that the noise from the atmosphere is dominant. 
The sensor index ranges between 1 and 76, corresponding to the wavefront Zernike coefficients 
Z4-Z22 at the four corners. 
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